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Abstract —We design a new secure transmission scheme in 
the relay wiretap channel where a source communicates with a 
destination through a decode-and-forward relay in the presence 
of spatially random-distrihuted eavesdroppers. For the sake of 
practicality, we consider a general antenna configuration in which 
the source, relay, destination, and eavesdroppers are equipped 
with multiple antennas. In order to confuse the eavesdroppers, we 
assume that both the source and the relay transmit artificial noise 
signals in addition to information signals. We first derive a closed- 
form expression for the transmission outage probability and an 
easy-to-compute expression for the secrecy outage probability. 
Notably, these expressions are valid for an arbitrary number of 
antennas at the source, relay, and destination. We then derive 
simple yet valuable expressions for the asymptotic transmission 
outage probability and the asymptotic secrecy outage probability, 
which reveal the secrecy performance when the number of anten¬ 
nas at the source grows sufficiently large. Using our expressions, 
we quantify a practical performance metric, namely the secrecy 
throughput, under a secrecy outage probability constraint. We 
further determine the system and channel parameters that 
maximize the secrecy throughput, leading to analytical security 
solutions suitable for real-world deployment. 

Index Terms —Physical layer security, wiretap channel, relay, 
secrecy outage, stochastic geometry, artificial noise. 


I. Introduction 

S ECURITY is a vital issue in wireless communication 
networks since data transmissions over the shared physical 
medium are inherently vulnerable to potential eavesdropping. 
Traditionally, security in wireless communication networks 
is realized by cryptographic techniques applied to the upper 
layers utilizing secret keys. The secrecy provided by such 
techniques is achieved under the assumption of hnite compu¬ 
tational capability at the eavesdroppers. However, this assump¬ 
tion cannot be easily satisfied with the rapid and continuous 
growth of the computational capability of modern proces¬ 
sors, which makes the traditional cryptographic techniques 
increasingly weak. Moreover, the ever-expanding size of de¬ 
centralized wireless networks introduces significant challenges 
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to key distribution and management. Against this backdrop, 
physical layer security has been proposed as a complementary 
technique to traditional cryptography, due to its benehts in 
enhancing the secrecy level of wireless communications by 
direct exploiting the randomness offered by wireless channels 
[1,2]. In seminal studies, e.g., [3], it was established in a 
single-input single-output wiretap channel that secrecy can 
only exist when the wiretap channel between the source and 
the eavesdropper is a degraded version of the main channel 
between the source and the legitimate receiver. This result was 
later generalized to the case where the main channel and the 
wiretap channel are independent [4]. 

Deploying multiple antennas at the source and/or the legiti¬ 
mate receiver has been shown to effectively boost the physical 
layer security of wiretap channels [5-17]. The effectiveness of 
multiple antennas relies on the use of secure multi-input multi¬ 
output (MIMO) techniques, such as beamforming [5-9], arti¬ 
ficial noise (AN) [10-13], and transmit antenna selection [14- 
17]. In the MIMO setting, the presence of randomly distributed 
eavesdroppers has been recently investigated [18-22]. In order 
to statistically characterize the secrecy performance of such 
scenarios, stochastic geometry and random geometric graphs 
are often used to model the locations of spatially random- 
distributed nodes. With such modeling, [18] investigated the 
throughput of large-scale decentralized wireless networks with 
physical layer security constraints. Considering the path loss 
as the sole factor affecting the received signal-to-noise ra¬ 
tios (SNRs) at the legitimate receiver and the eavesdropper, 
[19] examined the secrecy rate in cellular networks. In [20] 
and [21], the secrecy rate achieved by linear precoding was 
analyzed for the broadcast channel and the cellular network, 
respectively. In [22], the impact of AN was investigated. 

The above works [5-22] examine physical layer security in 
point-to-point MIMO systems. Cooperative relaying, on the 
other hand, is another promising and widely-adopted tech¬ 
nique that efficiently improves the coverage and reliability of 
wireless networks [23,24]. In order to enhance physical layer 
security in relay wiretap channels, a variety of approaches 
have been investigated such as cooperative beamforming [25- 
28], relay selection [29,30], and cooperative jamming [31, 
32]. However, a common limitation of [25-32] is that they 
only considered fixed locations of eavesdroppers. This leaves 
open the problem of designing relay-aided secure transmission 
schemes for the scenario where the locations of eavesdroppers 
are spatially randomly distributed. 
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In this work we design a new relay-aided secure trans¬ 
mission for the relay wiretap channel. In such a channel, 
the communication between the source and the destination 
is aided by a decode-and-forward (DF) relay and overheard 
by multiple spatially random-distributed eavesdroppers. We 
focus on the general scenario where the source, the relay, the 
destination, and the eavesdroppers are equipped with multiple 
antennas, which stands as a major advancement over the 
previous studies on securing the relay wiretap channel [25- 
32]. In order to confuse the eavesdroppers, we assume that in 
the secure transmission the source and the relay transmit AN 
signals together with information signals in the first hop and 
the second hop, respectively'. The contributions made by this 
work are summarized as follows: 

1) We derive a closed-form expression for the transmission 
outage probability and an easy-to-compute expression 
for the secrecy outage probability. Notably, both expres¬ 
sions are independent of realizations of channels and 
valid for an arbitrary number of antennas at the source, 
relay, and destination. Moreover, these expressions serve 
as the key results that enable us to explicitly characterize 
the secrecy throughput of the considered relay wiretap 
channels. 

2) We derive simple yet valuable expressions for the 
asymptotic transmission outage probability and the 
asymptotic secrecy outage probability. These expres¬ 
sions quantify the secrecy performance in the regime 
where the number of antennas at the source becomes 
sufficiently large. Based on our analysis, we find that 
the asymptotic transmission outage probability is de¬ 
termined by the average SNR of the relay-destination 
channel only. We also find that the asymptotic secrecy 
outage probability approaches a certain value which is 
independent of the number of antennas at the source. 

3) We determine the transmission parameters, i.e., the 
wiretap code rates and the power allocation factors, that 
maximize the secrecy throughput of the considered relay 
wiretap channels under a secrecy outage probability 
constraint. Moreover, we demonstrate the effectiveness 
of the determined transmission parameters on maximiz¬ 
ing the secrecy throughput. Furthermore, we evaluate 
the impact of the system parameters, e.g., the number 
of antennas and the density of eavesdroppers, on the 
secrecy throughput. 

Beyond the above contributions, we provide some pivotal 
insights into the practical design of secure transmission. First, 
we show that the AN signals from the source play a more dom¬ 
inant role in securing the transmission in the considered relay 
wiretap channels than the AN signals from the relay. Second, 
we show that adding extra antennas at the source significantly 
increases the maximum secrecy throughput, but does not 
decrease the secrecy outage probability always. Third, we find 
that in order to achieve the maximum secrecy throughput, the 

*An initial study of a much simpler system model is given in [33] where 
the relay, the destination, and the eavesdroppers are all equipped with a single 
antenna and AN signals are transmitted by the source in the first hop only. 
This simplified system configuration allowed for analytical tractability at the 
expense of significant sub-optimality. 




Fig. 1. Illustration of a relay wiretap channel in the presence of spatially 
random multi-antenna eavesdroppers. 


source needs to allocate a higher power to AN signals whereas 
the relay needs to allocate a lower power to AN signals 
when the antenna number at the source increases. Fourth, we 
find that the maximum secrecy throughput increases when the 
eavesdroppers are more dispersed. 

The rest of the paper is organized as follows. Section II 
describes the relay wiretap channel considered in the paper. In 
Section III, we derive expressions for the outage probabilities 
of the considered relay wiretap channel. The characterization 
and maximization of the secrecy throughput are also provided 
in Section III. Numerical results and related discussions are 
presented in Section IV. Finally, Section V draws conclusions. 

Notations: Column vectors (matrices) are denoted by bold¬ 
face lower (upper) case letters. Conjugate transpose is denoted 
by ■ The determinant of a matrix is denoted by det (•). 
Complex Gaussian distribution is denoted by CM. A zero 
matrix and an identity matrix of appropriate dimension are 
denoted by 0 and I, respectively. Statistical expectation is 
denoted by E. The Frobenius norm of a vector or a matrix 
is denoted by || • ||. 

II. Multi-Antenna Relay Wiretap Channel 

We consider a relay wiretap channel, as depicted in Fig. 1, 
where a source (S) communicates with a destination (D) with 
the aid of a relay (R) in the presence of multiple spatially 
random eavesdroppers. In this channel, the source, the relay, 
the destination, and each eavesdropper are equipped with Ng, 
Nr, Nfi, and N^. antennas, respectively. We denote as 
the Nr X Ns channel matrix from the source to the relay and 
denote as the Nd x Nr channel matrix from the relay to 
the destination. We consider that all the channels are subject to 
independent and identically distributed (i.i.d) Rayleigh fading. 
We also consider a quasi-static block fading environment in 
which all the channel coefficients remain the same within one 
time slot. We assume that the channel state information (CSI) 
between the source and the relay and the CSI between the 
relay and the destination are known at the source, while the 
CSI from the eavesdroppers is not known. We also assume 
that Ns > Ne, mimicking the case where the source is a 
base station (BS) with a large number of antennas, while 
the eavesdroppers are mobile users with a limited number of 
antennas. We further assume that the destination is located 
remotely away from the source such that the destination 
cannot receive signals from the source directly. All the nodes 
operate in a half-duplex mode such that each node cannot 
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transmit and receive simultaneously. We denote dsr and drd 
as the source-relay distance and the relay-destination distance, 
respectively, and denote 77 as the path loss exponent. The 
locations of the eavesdroppers are modeled as a homogeneous 
Poisson Point Process (PPP) $ with density A [18-22], which 
represents the case where the eavesdroppers are mobile users 
in a decentralized network [34]. We clarify that the source, 
the relay, and the destination do not belong to $. 


A. Transmission of Artificial Noise Signals 

We now detail the transmission scheme between the source 
and the destination. In this scheme we assume that both the 
source and the relay transmit AN signals together with the 
information signals. This scheme utilizes two time slots. In 
the first time slot, the source transmits information signals 
and AN signals to the relay, referred to as the first hop 
transmission. We assume that the relay adopts maximum-ratio 
combining (MRC) [35-37] to process the received signals 
in order to maximize the received SNR. In the second time 
slot, the DF relay transmits the re-encoded signals and AN 
signals to the destination, referred to as the second hop 
transmission. We assume that the destination also adopts MRC 
to process the received signals. In the second time slot, it 
is assumed that the source transmits AN signals to further 
confuse the eavesdropper. We clarify that both the first hop 
transmission and the second hop transmission are overheard 
by the eavesdroppers. 

In the first hop transmission, the signal transmitted by the 
source is given by 

xs = Witi, (1) 


where Wi denotes the Ng x Ng beamforming matrix at the 
source and ti denotes the combination of the information 
signal and the AN signal at the source. To transmit xs, we 
first design Wi as 

Wi = [ws Wsan] , (2) 


where ws is used to transmit the information signal at the 
source and Wsan is used to transmit the AN signal at the 
source. The aim of Wi is to degrade the quality of the 
received signals at the eavesdroppers. By transmitting AN 
signals through Wi, together with the fact that the relay adopts 
MRC to process the received signals from the source, we 
ensure that the quality of the received signals at the relay 
is free from AN interference. In designing Wi, we choose 
Ws as the eigenvector corresponding to the largest non-zero 
eigenvalue of H^Hgr, denoted by We then choose 

Wsan as the remaining Ng — 1 eigenvectors of H^Hgr. Such 
design ensures that W 1 is a unitary matrix. We then design t 
as 


ti — 


h 

tsAN ’ 


(3) 


where ts denotes the information signal at the source and tsAN 
is an (Ng — 1) x 1 vector of the AN signal at the source. 
We define fig, 0 < fig < 1, as the fraction of the power 
allocated to the information signal at the source. As such, we 


have E [|fsP] = fig and E [tsANt|^N] = Based 

on (1), (2), and (3), the received signal at the relay in the first 
hop transmission is expressed as 

yr = \J Pgd7fiti.gr (wsts + WsANtsAN) + n^, (4) 

where Pg denotes the transmit power at the source and 11 ^ 
denotes the thermal noise at the relay, the elements of which 
are assumed to be i.i.d complex Gaussian random variables 
with zero mean and variance a^, i.e., ~ CAf cr^Ijv,,). 

We note that AN signals in (4) can be canceled at the relay 
by applying MRC. 

We next express the received signal at a typical eavesdrop¬ 
per located at i, z € 4>, in the first hop transmission as 

= \JPi>dfl''iigi (wsfs + WsANtsAN) + Wi, (5) 

where Hgi denotes the Ne x Ng channel matrix from the 
source to the typical eavesdropper located at z, dgi denotes 
the distance between the source and the typical eavesdropper 
located at z, and 11^1 denotes the thermal noise vector at the 
typical eavesdropper located at z, the elements of which are 
assumed to be i.i.d complex Gaussian random variables with 
zero mean and variance i.e., nn ~ CM (O^Ve, cr^ilAZe)- 
In the second time slot, the DF relay first decodes the 
received signals from the source. If the received signals are 
successfully decoded, the relay retransmits the re-encoded sig¬ 
nals and AN signals to the destination. The signals transmitted 
by the relay is given by 

Xr = W2t2, (6) 


where W 2 denotes the Nr x Nr beamforming matrix at the 
relay and t 2 denotes the combination of the information signal 
and the AN signal at the relay. Similar to Wi and ti, we 
design W 2 and t 2 as 


W 2 = Wr Wran , (7) 


and 


t 2 


tRAN ’ 


( 8 ) 


respectively. In (7), wr is used to transmit the information 
signal at the relay and Wran is used to transmit the AN 
signal at the relay. In designing W 2 , we choose wr as 
the eigenvector corresponding to the largest eigenvalue of 
denoted by We then choose Wran as the 

remaining A^^ —1 eigenvectors of This design ensures 

that the quality of the received signals at the destination is free 
from AN interference when the destination applies MRC to 
process the received signals. In ( 8 ), Ir denotes the information 
signal at the relay and tRAN is an {Nr — 1 ) x 1 vector of 
the AN signals at the relay. We define fir, 0 < fir < 1, 
as the fraction of the power allocated to the information 
signals at the relay. As such, we have E [|fRp] = fir and 
E [tRANtRAN] = T^'^Nr-i- According to ( 6 ), (7), and ( 8 ), 
we express the received signal at the destination in the second 
hop transmission as 


yd 



-n 

rd 


Hrd (wr^r -f WRANtRAx) + Hd, 


(9) 
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where Pr denotes the transmit power at the relay and 
denotes the thermal noise at the destination, the elements of 
which are assumed to be i.i.d complex random variables with 
zero mean and variance i.e., ~ CjV 

We note that AN signals in (9) can also be canceled at the 
destination by applying MRC. 

In order to further confuse the eavesdroppers in the second 
hop transmission, we assume that the source transmits AN 
signals using transmit power Pg. We denote the AN signals 
from the source in the second hop transmission as xan, 
the elements of which follow the i.i.d zero mean complex 
Gaussian distribution. We assume that xan has unit power such 
that E [xan^an] = '^nJNs. We next express the received 
signal in the second hop transmission at a typical eavesdropper 
located at i, i G $, as 

yf'’ =\JPrd~^'B.ri (wRfR + WRANtRAx) 

+ \J Psd“’'HsiXAN + ni2, ( 10 ) 


respectively, where 

= ^?^^.rf.7H,.WsANWs"ANHf! + (14) 

and 

=^^^.7"H„WRANW^ANH^i 

+ + (15) 

We assume that the eavesdroppers are non-colluding, indicat¬ 
ing that each eavesdropper decodes her own received signals 
from the source and the relay without cooperating with other 
eavesdroppers. We also assume that the source and the relay 
use different codebooks. As such, the transmitted signals from 
the source and the transmitted signals from the relay cannot 
be jointly processed at each eavesdropper. Based on (12) and 
(13), we express the equivalent SNR at the eavesdroppers as 

Fa; = max{max{7si,7ri}} . (16) 


where denotes the x Nr channel matrix from the 
relay to the typical eavesdropper located at i, dri denotes 
the distance between the relay and the typical eavesdropper 
located at i, and 11^2 denotes the thermal noise vector at a 
typical eavesdropper located at i, the elements of which are 
assumed to be i.i.d complex Gaussian random variables with 
zero mean and variance ctIj, i.e., ni 2 ~ CAf (0N^,af2lNe)- 

B. Formulation of Received Signal-to-Noise Ratios 

We first focus on the equivalent instantaneous SNR at the 
destination. Recall that both the relay and the destination 
apply MRC to process received signals. We express the MRC 
combiner at the relay in the first hop transmission as = 
TiTT—and express the MRC combiner at the destination in 

the second hop transmission as = ||h^ • Using and 

Vd, we express the instantaneous SNR at the relay in the first 
hop transmission and the instantaneous SNR at the destination 
in the second hop transmission as 73 ^ = 

'Ird = respectively. As per the rules of the DF 

protocol, we express the equivalent end-to-end SNR from the 
source to the destination as [23] 

r_D = min{7sr,7rd} • (11) 

We now focus on the equivalent SNR at the eavesdroppers. 
In order to maximize the probability of successful eavesdrop¬ 
ping, we assume that the eavesdropper utilizes the minimum 
mean square error (MMSE) combining to process the received 
signals within two time slots. As per the rules of the MMSE 
combining, we express the instantaneous SNR at a typical 
eavesdropper located at i in the first hop transmission and 
the second hop transmission as 

H«K7H,,ws, (12) 


III. Secrecy Performance Analysis 
In this section, we analyze the secrecy performance achieved 
by the transmission scheme detailed in Section II. We first 
derive a closed-form expression for the transmission outage 
probability and an easy-to-compute expression for the secrecy 
outage probability, both of which are valid for an arbitrary 
number of antennas at the source, relay, and destination. We 
then derive simple yet valuable expressions for the asymptotic 
transmission outage probability and the asymptotic secrecy 
outage probability, both of which are valid for a sufficiently 
large number of antennas at the source, i.e., Ng —> 00 . We 
further describe in detail how the secrecy throughput of the 
relay wiretap channel is quantified and how the maximum 
secrecy throughput is obtained under a secrecy outage proba¬ 
bility constraint. 


A. Preliminaries 


In this subsection, we present the statistics of ■jgr, 'jrd, Isi, 
and 7 j.j, which will be used to derive the outage probabilities. 
We first focus on the cumulative distribution functions (CDEs) 
of 7 sj. and 7 ^^. To this end, we introduce several new notations 
as follows: ui = min(A^s, A^^.), vi = ma.x{Ns, Nr), ti = 
vi — ui, U 2 = min (Nr, Nf), V 2 = max{Nr, Nd), and t 2 = 
V 2 — U 2 - We then obtain the CDE of 7 sr as [38] 


Fj.r ( 7 ) 


r„i (ui)r„^ (ui) ’ 


(17) 


where 


5 ( „ Z ) is a Ml X ui matrix with (i,j)th entry, 

\ Hsisr J 

(a)’ 7 


In (18), 7 (•) denotes the incomplete gamma function, defined 
as [39, Eq. (8.352)] 


fc-i 


1 - exp (-x) ^ — 


2=0 


and 


7„ = /Sr^rfi^WR 


(13) 


7 (fc,x) = r(A:) 


(19) 
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for integer k, where F (•) denotes the gamma function, defined 
as r {k) = [k — 1)! for integer k [39, Eq. (8.339)], gi {i,j) = 
ti + i + j — Isr = Psdj^<7~'^, and 

n 

Tm (n) = n r (m - t + 1 ). ( 20 ) 

2=1 


Similarly, we obtain the CDF of ^rd as 


iu2) (W 2 ) ’ 


where 0 Z j is a W 2 x ^2 matrix with (i,j)th entry, 


of the wiretap code, and R^, denotes the redundancy rate of the 
wiretap code revealing the cost of preventing eavesdropping. 
We also assume that the source and the relay use the same 
{Rb^ Re) to transmit, but with different codebooks. As such, 
we define that the transmission outage event occurs when 
Cb < Rb- In this event, the received signals at the destination 
are not reliably decoded. We also define that the secrecy 
outage event occurs when Ce > Re- In this event, the 
eavesdropper is able to decode the transmitted signals and 
secrecy is compromised. 

Based on the definition of the transmission outage event, we 
define the transmission outage probability as the probability 
that the equivalent instantaneous SNR at the destination is 
less than Tb = 2^'’ — 1. Mathematically, Pto is formulated as 


% () =7(52(*,j),-T^ ) , ( 22 ) 

KPrlrdJ \ PrlrdJ 

where 52 (*, j) = ^2 + * + J - 1 , and 7 ^^ = Prd~^a^'^. 

With the aid of [40], we express the CDF of jsi as 

F 1 1 ( ^ Y~" 

X "E TZ" ') (-17)’ (23) 

where 7 ^^ = Psd~l'cT--^ and m = ^ ^ and express 

the CDF of jri as 

Fjri ( 7 ) 

(1 + K27)^’'"^ (1 + K37)^'’ P ("l) 


Nr - 1 


(K 27 )” E 


(«37)*, (24) 


where 7 ^^ = Prdrl'cr^^^, K 2 = wkS 


l3r(Nr-l) ’ -3 jdrPrNsd'l 


Pto = Pr (F^ < Tb). (27) 

Using (11), (17), and (21), we re-express the transmission 
outage probability in (27) as 

Pto = Pr (min { 7 ^^, ^rd} < u) 

= 1 - (1 - {Tb)) (1 - {Tb)) 

r«i (ui) (Wl) (M 2 ) (M 2 ) 

_ •^°>(-(ia:))'^°>('^(ftt 7 )) 

(mi) (mi) ^122 (M 2 ) ^172 (M 2 ) ■ 

Based on the definition of the secrecy outage event, we 
define the secrecy outage probability as the probability that F e 
is larger than Te = 2^® — 1. Mathematically, Pgo is formulated 
as 

Pso = Vr{TE>Te). (29) 

According to (16), (23), and (24), we derive an easy-to- 
compute expression for the secrecy outage probability in the 
following theorem. 

Theorem 1: The secrecy outage probability of the relay 
wiretap channel is derived as 


B. Outage Probabilities 

In this subsection, we define the transmission outage event 
and the secrecy outage event and then characterize their 
probabilities. We first denote Cb as the instantaneous capacity 
between the source and the destination. According to (11), Cb 
is given by 

Cb = kog^{l + TD), (25) 

where the presence of the factor 1/2 is due to the fact that 
two time slots are used in the transmission. We also denote 
Ce as the instantaneous capacity between the source and the 
eavesdropper. According to (16), Ce is given by 

Ce = ^ log 2 (1 + Te) . (26) 

We assume that the wiretap code is adopted in the trans¬ 
mission. We denote {Rb,Re) as the parameter pair for the 
adopted wiretap code, where Rb denotes the transmission rate 


P«o = l-exp(-2A(Ji+J2-(^3)), (30) 



N 2 and Ns are given by (32) and (33), respectively, shown at 
the top of the next page. In (32) and (33), we have ip {0) — 

^77 ksr + dl, - 2dsrdsi cos 9) ^ . 

Proof: See Appendix A. ■ 

We find that Theorem 1 provides an easy-to-compute tool 
for efficiently evaluating the secrecy outage probability. Al¬ 
though J 2 and J-i for general rj cannot be obtained in closed- 
form, they can be easily calculated since only a double integral 
is involved in J 2 and 
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n TT 

d, 


exp(-V> (6*)) 
Y 1 _L 


Ne 


1 (m) 

m—l ^ ' 


Ne — m 

E 

n—O 


Nr - 1 


Ne — Tn—n 


{K2TeT E 


i=0 


Ns\ ( Psi’ie) 


ddsidd, 


(32) 


Js = (1 + KlTe) (1 + K2Te) 




^ exp(-V> (6*))^^ ^ 1 ^E * ^ 


V 


N, 


n—0 


Ne —m—r 


9=0 


{n2TeT E 


i=0 


Ns\ ( Psi’iO) V 


dd^jdO. 


(33) 


C. Asymptotic Outage Probabilities 

In this subsection, we examine the asymptotic behavior of 
the outage probabilities as Ng oo. The obtained asymptotic 
results are particular valuable for large-scale MIMO systems 
where the source (or equivalently, the BS) is equipped with 
a sufficiently large number of antennas. We first present the 
expression for the asymptotic transmission outage probability 
in the following corollary. 

Corollary 1: The asymptotic transmission outage probabil¬ 
ity when Ng ^ oo is given by 


det 


poo _ 

^to 




{U2}r„, (U2)' 

Proof: We express the asymptotic transmission outage 
probability when Ng ^ oo S.S 


Pff= lim 

tV.->oo 


Pf. 


We note that 


lim 




= 0 . 


(35) 


(36) 


n7Nco r„^ (ui) (mi) 

Substituting (36) into (35) yields the result. ■ 

According to Corollary 1, we find that the asymptotic 
transmission outage probability is solely determined by 
when Ng oo. This finding is due to the fact that jgr oo 
when Ng —)■ oo. As such, we conclude that the probability 
that Td is less than Tb when > oo is determined by the 
link quality of the relay-destination channel only. 

We next present the asymptotic secrecy outage probability 
when oo in the following corollary. 

Corollary 2: The asymptotic secrecy outage probability 
when Ng ^ oo is given by 


P“ = 1 - exp (-2A (Jr + Jr - J 3 “)), 


(37) 


where 


1 \TeafJ V Ps 


r (p) ^ 


P=1 


9=0 


Tiq+1) ’ 


(38) 


J 2 °° and J'f° are given by (39) and (40), respectively, shown 
at the top of the next page. 

Proof: We express the asymptotic secrecy outage proba¬ 
bility when —>■ oo as 


poo 

^ SO 


lim Pgo. 

Ng—^OO 


(41) 


We note that 


lim f 1 + 




j3g{Ng — l) 

and 


= exp - 


1-r 


lim 

Ng —^OO 


We also note 


lim ( 1 + 

Ng—foo 


(iM 

PgiNg-l)J Tiq+1) ■ 


Ng-1\ f (l-r)re 


Psf jO) 

NsCpg^crh 


-N, 


= exp 


and 


Ng\ f Pgf (9) 


tvl™oo \ I J \ Ngdlal 


Psf jO) 

dlcr^2 


( PNi+ y 

V J 

r(( + i) ■ 


(42) 


(43) 


(44) 


(45) 


Substituting (42), (43), (44), and (45) into (41) yields (37), 
which completes the proof. ■ 

According to Corollary 2, we find that the asymptotic 
secrecy outage probability approaches a certain value that is 
independent of Ng when Ng —^ oo. This reveals that adding 
extra transmit antennas at the source does not always decrease 
the secrecy outage probability. 


D. Secrecy Throughput 

So far, we have derived the exact and the asymptotic 
transmission outage probability and secrecy outage probability 
of the considered relay wiretap channel. Our derived expres¬ 
sions are valid for given Rb, Re, fig, and fir- A question 
then naturally arises: “How do we determine the optimal 
{Rff, Rff, j3f°, fi*°) that achieves the maximum secrecy per¬ 
formance of this relay wiretap channel, under a secrecy 
outage probability constraint?” Our answer to this question 
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(40) 


demonstrates the usefulness of our analytical expressions in a 
wider sense. It shows how the expressions can be embedded 
and utilized in a complex optimization problem, presenting 
a solution that can be determined in a much faster manner 
than would otherwise be possible. A concrete example of 
this usefulness in an operational sense, would be the dynamic 
and real-time determination of the optimal system parameter 
settings for a given secrecy performance metric. 

To answer the question, we utilize our derived expressions 
directly to characterize the secrecy performance of the relay 
wiretap channel. As a specific example, we consider the metric 
termed the secrecy throughput introduced by [10]. This perfor¬ 
mance metric quantifies the average confidential information 
rate when the source transmits. The secrecy throughput for the 


relay wiretap channel is given by 

Ts = {Rb - Re) (1 - Pto) ■ (46) 

The maximization problem is accordingly formulated as 

max Ts, (47 a) 

Rb,Re,0s:/3r 

s.t. Pso < y, (47b) 

0 <Re < Rb, (47c) 

0 </3s < 1,0 </3^ < 1. (47d) 


In the following, we describe in detail how the maximum 
secrecy throughput is obtained by judiciously selecting the 
transmission parameters^. To this end, we solve the maxi¬ 
mization problem in (47) in two steps. First, we fix power 
allocation factors Ps and Pr, and choose the wiretap code 
rates pair, that maximizes the secrecy throughput. 

Accordingly, the maximum secrecy throughput achieved by 
{Rl,Rl) for given Ps and Pr is defined as T*. Second, we 
choose the wiretap code rates as well as the power allocation 
factor, {Rl°, Rl°, P*°, P*°), that jointly maximizes Ts. The 
details of these two steps are presented as follows: 

1) {Rl,R*) for given Ps and Pr-' The wiretap code rates 
pair, {Rl,Rf), that maximizes Tg for given Ps and Pr is 

^We note that the maximum secrecy throughput is formally a local max- 
imum, since we numerically search the transmission parameters that jointly 
maximize the secrecy throughput only within anticipated ranges. 


determined as 

{Rl,Rl) = argmax Ts, (48a) 

s.t. Pso < p, (48b) 

0 < Re P: Rb- (48c) 

Taking the first-order derivative of P^o with respect to Re, 
we confirm that dPso/dRe < 0, which indicates that Pso 
monotonically decreases as Re increases. As such, the value 
of P* satisfying (48b) is the value of P* that satisfies the 
secrecy outage probability constraint, i.e., Pgo (Rl) = p. 
We then confirm that dPto/dRb > 0, which shows that Pto 
monotonically increases as Re increases. As such, we note that 
Ts —)■ 0 as Rb increases such that Pto —> 1. Defining as 
the value of Rb that satisfies Pto (Rb) = 1, we rewrite (48) as 

{Rl,R*e) = argmax T,, (49a) 

s.t. Pso < p, (49b) 

R*e<Rb< (49c) 

Although a closed-form solution for (P^,P*) is mathemati¬ 
cally intractable, we are able to find the values of (P ^, Rl) in 
a numerical way. 

2) {Rl°, P*°, P*s°,P*°): The wiretap code rates and power 
allocation factors which jointly maximizes T, in (46), 
{Rl°, Rl°, P*°, P*°), is determined as 

{Rl°,RT,p*s°,P*°) = argmax Ts, (50a) 

s.t. Pso < p, (50b) 

0 < Pe < Rb, (50c) 

0 </3s < 1,0 << 1. (50d) 

Using (28) and (30), we are able to solve (50) numer¬ 
ically. Specifically, we first select the value of (P^,P*) 
satisfying (49) for each value of Pa and Pr. This leads to 
the secrecy throughput with (P^,P*), denoted by T* = 
(Pj — Pg) (1 — Pto). We then select the value of P*° and 
the value of Pr° that maximize T* for 0 < /3s < 1 and 
0 < /3r. < 1. Accordingly, the value of {RfRf) associated 
with Pl° and /3*° is defined as (P^°, Pe°). Finally, the maxi¬ 
mum secrecy throughput achieved by {Rl°, Rl°, P*°, P*°) is 
defined as T*°. 
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Fig. 2. Pto versus for different values of 7 ^, with 77 = 4, Ns = 4, Nr = 2, 
Nfi = 2, Ne = 2, ps = pr = 0.5, A = 0.01, and dsr = drd = 10. 



T. (dB) 

Fig. 3. Pso versus Te for different values of 7 g with r; = 4, Ns = 4, Nr = 
2, Nd = 2, Ne = 2, 13s = /3r = 0.5, A = 0.01, and dsr = drd = 10. 


We note that our expressions for the outage probabilities 
can also be implemented in other performance metrics. For 
example, we can utilize our expressions to characterize the 
average secrecy rate of the relay wiretap channels in the 
presence of spatially random eavesdroppers. 

IV. Numerical Results 

In this section, we present numerical results to validate 
our analysis of the outage probabilities. We also examine 
the impact of transmission parameters (e.g., Rj,, fig, Pr) and 
system parameters (e.g., Ng and A) on the secrecy throughput 
of the relay wiretap channel. Throughout this section we 
concentrate on the practical example of a highly shadowed 
urban area with p = 4. 

We hrst demonstrate the accuracy of the transmission outage 
probability and the secrecy outage probability using Monte 
Carlo simulations. In Fig. 2, we plot Pto versus rj, for different 
values of 7 f, with Ng = 4, = 2, Nd = 2, Ng = 2, 

Ps = Pr = 0.5, A = 0.01, and dgr = drd = 10. In this 
figure, we consider 7 ,,^ = j^d = 76 - We first see that the 
analytical curves, generated from (28), precisely match the 
simulation points marked by red circles, which demonstrates 
the correctness of our expression for Pto in (28). Second, 
we see that Pto increases monotonically as Tb increases for a 
given 7 f,. This reveals that the transmission outage probability 
increases when the transmission rate of the wiretap code 
increases. We further see that Pto decreases as 7 ^, increases for 
a given Tb- This reveals that the transmission outage probability 
reduces when the source and the relay use a higher power to 
transmit for a fixed tj,. 

In Fig. 3, we plot Pgo versus Tg for different values of 7 ^ 
with Ng = A, Nr = 2, Nd = 2, Ng = 2, and jig = j3r = 
0.5. In this hgure we consider ^gr^^rl^'ii = lrdP\l^'i 2 = 
7 g. We see an excellent match between the analytical curves 


generated from (30) and the simulation points marked by red 
circles, demonstrating the correctness of our expression for Pgo 
in (30). We then see that Pgo decreases monotonically as Tg 
increases for a given 7 ^, which shows that the secrecy outage 
probability decreases when the redundancy rate of the wiretap 
code increases. We further observe that Pgo increases as 7 g 
increases. This is due to the fact the eavesdroppers receive 
signals from both the source and the relay. It follows that 
increasing the transmit power at the source and the relay leads 
to an improved received SNR at the eavesdroppers. 

In the following, we examine the impact of the transmission 
parameters on the secrecy throughput that is characterized by 
the derived outage probabilities. We first examine the impact of 
Rb on Tg. In Fig. 4, we plot Tg versus Rb for different values 
of Ng with i?g and fixed Nd, Ng, jdg, and j3r- We hrst observe 
that there exists a unique R^ that maximizes Tg for given fdg 
and Pr- We also observe that the maximum Tg for given Pg and 
Pr, i.e., Tg, increases as Ng increases. This shows that adding 
extra transmit antennas at the source signihcantly enhances the 
secrecy performance of the relay wiretap channel. 

We next examine the impact of Ps and Pr on T*. In Fig. 
5, we plot T* versus Pg for different values of Ng with a 
hxed pr- For each point of T*, we choose {Rl,R*) that 
maximizes Tg for the corresponding Pg. We hrst observe 
that there exists a unique P*° that maximizes T*. We then 
observe that the maximum T* for a hxed Pr increases as 
Ng increases. Furthermore, we observe that the value of 
Pg° slightly decreases as Ng increases, which shows that in 
order to maintain the maximum secrecy throughput, the power 
allocated to AN signals at the source needs to be increased as 
the number of antennas at the source increases. 

In Fig. 6 , we plot T* versus pr for different values of Ng 
with a hxed Pg. Similar as Fig. 5, for each point of T*, we 
choose (i?^,i?g) that maximizes Tg for the corresponding Pr. 
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Fig. 4. Tg versus Ri, for different values of Ng with rj = A, = 2, 
Na = 2, Ng = 2, 0g= pr = 0.5, A = 0.01, dgr = = 10, tp = 0.4, 

7 f, = 10 dB, and 7i,/7e = 20. 



Fig. 5. Tg versus /3s for different values of Ng with rj = A, Nr = 2, 
Ni = 2, Ng = 2, pr = 0.5, A = 0.01, dgr = drd = 10, p = 0.4, 
7 t, = 10 dB, and 7 t,/ 7 e = 20. 

First, we see a unique /3*° that maximizes T*. Second, we see 
that the maximum T* for a fixed fig increases as Ng increases. 
Additionally, we note that the value of /3*° increases as Ng 
increases, demonstrating that a lower power is needed to be 
allocated to AN signals at the relay in order to achieve the 
maximum secrecy throughput when the antenna number at the 
source increases. 

Finally, we examine the impact of Ng and A on T°*. In Fig. 
7, we plot T*° versus Ng for different values of A. For each 
point of T*°, we choose , Pr, Rl°, Rl°) that maximizes 



Fig. 6. Tg versus 0r for different values of Ng with ij = A, Nr = 2, 
Nd = 2, Ng = 2, Pg = 0.5, A = 0.01, dgr = drd = 10, ip = 0.4, 
7 f, = 10 dB, and 7(,/7e = 20. 



Ns 


Fig. 7. Tg° versus Ng for different values of A with rj = A, Nr = 2, Nd = 
2, Ng = 2, dgr = drd = 10, tp = 0.4, 7(, = 10 dB, and 7i,/7e = 20. 


Tg. We first observe that T/° increases as Ng increases. This 
observation is consistent with the observation in Fig. 4. We 
then observe that T/° decreases as A increases. This is due 
to the fact that more eavesdroppers exist as A increases. The 
increasing number of eavesdroppers increases the value of i?| 
that satisfies the secrecy constraint. 

V. Conclusion 

We designed the secure transmission that maximizes the 
secrecy throughput of the generalized relay wiretap channel in 
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the presence of spatially random multi-antenna eavesdroppers. 
In the transmission we assumed that both the source and 
the relay transmit AN signals with information signals in 
order to confuse the eavesdroppers. Considering the use of 
the decode-and-forward relaying protocol, we first derived a 
closed-form expression for the transmission outage probability 
and an easy-to-compute expression for the secrecy outage 
probability. We then derived simple yet valuable expressions 
for the asymptotic transmission outage probability and the 
asymptotic secrecy outage probability when the number of 
antennas at the source becomes sufficiently large. Using our 
derived expressions, we characterized the secrecy throughput 
of the the relay wiretap channel and then determined the 
transmission and system parameters that achieve the maximum 
secrecy throughput. Finally, we evaluated the impact of these 
parameters on the secrecy throughput. 


Jo Jo (1 -I- KlTe)^‘~^ 

We , / 2 \ m-1 Ne-m e,, , 

/'v 7 "- 


I J \/3rPrNsdl 


Appendix A 
Proof of Theorem 1 

According to (16), (23), and (24), we re-express (29) as 
Pso = 1 - Pr{r£; < Te} 


= 1 - Pr max {max { 7 ^*, jri}} < A 
ie$ 


= 1-E$ Pr{max{ 7 si, 7 H} < Te} 


= 1 - E$ (Te) (Te) 


1 - exp (-2A {Ji +J 2 - Jo)). 


(•00 eTT / exp 

/ dsi I 




(\ We / 2 \ P-1 

, _ \N,-1 r I R. p. 


N, - 1 


(KiTe)® ddsidO, (49) 




[1 + ^ (1 + 


P3 Te d^i 
drPrN.d'O 


We / 2 \ m — l We —m / -, , 

f PsredJ. .. 

£ (, 1 j \apmiJ 


In (48), the operation (a) can be justified by applying the 
probability generating functional (PGFL) for the PPP $, given 
by [41] 

]E<[. TT./(a;) =exp|-/' [1 - f (x)] Xdx\ , (52) 

Ue$ J I 4 r 2 j 


and by changing to polar coordinates. 

In order to proceed with our analysis we first derive Ji as 


^ 1 fN, - 1 


J, =7r (1 + NiTe) E E 


(NlTe)'^ 


We We-p 


AAidl dd. 
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where in (6) we use u = d^^, in (c) we use t = Yp 
and [d) follows from the definition of the gamma function. 
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Similarly, we derive J 2 as 
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where in (e) we apply the cosine formula. We further derive 

Js as 
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Substituting (53), (54), and (55) into (48), we obtain the 
desired result in (30), which completes the proof. 
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